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ABSTRACT 

Clusters of galaxies are revealing themselves as powerful sources of non 
thermal radiation in a wide range of wavelengths. In order to account for these 
multifrequency observations equipartition of cosmic rays (CRs) with the thermal 
gas in clusters of galaxies is often invoked. This condition might suggest a 
dynamical role played by cosmic rays in the virialization of these large scale 
structures and is now testable through gamma ray observations. We show here, 
in the specific case of the Coma and Virgo clusters, for which upper limits on 
the gamma ray emission exist, that equipartition implies gamma ray fluxes that 
are close or even in excess of the EGRET limit, depending on the adopted model 
of CR injection. We use this bound to limit the validity of the equipartition 
condition. We also show that, contrary to what claimed in previous calculations, 
the equipartition assumption implies gamma ray fluxes in the TeV range which 
can be detectable even by currently operating gamma ray observatories if the 
injection cosmic ray spectrum is flatter than E~ 2A . 

Subject headings: Galaxies:clusters:general — gamma rays:theory 
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1. Introduction 

Several non-thermal processes have recently been detected in clusters of galaxies from 
the extreme ultraviolet (EUV) radiation in excess of the thermal expectation to the soft 
X-rays detected by ROSAT and BeppoSAX and again to hard X-ray excesses and radio 
radiation. For the Coma cluster, by far the best studied cluster, a complete investigation of 
the soft excess can be found in (Lieu et al. 1996a, 1996b; Bowyer, Lampton & Lieu 1996; 
Fabian 1996; Mittaz, Lieu & Lockman 1998, Sarazin & Lieu 1998) while the detection of 
the hard excess above 20 keV is reported in (Fusco-Femiano et al. 1998). A recent review 
of the diffuse radio emission can be found in (Feretti et al. 1998) . Some clusters show only 
emission in some region of frequency and not in others. Also for this reason Coma gives the 
best possibility to make multiwavelengths studies. A review of the current status of the 
multifrequency observations of Coma and viable models for the non thermal radiation can 
be found in (Ensslin et al. 1998) . As stressed by Fusco-Femiano et al. (1998), if the hard 
X-ray excess is due to inverse compton scattering (ICS) off the photons of the microwave 
background, then the combined radio and hard X-ray observations of Coma imply a small 
value for the average intracluster magnetic field, of order B ~ 0.1 — 0.2fiG. Such a small 
value of the field requires large energy densities in electrons, and, as pointed out in (Lieu et 
al. 1999), CR energy densities comparable with the equipartition value are required. This 
conclusion is only weakly dependent on the specific model (primary or secondary) for the 
production of the electrons responsable for the radiation. In fact the need for large CR 
energy densities was recently confirmed by Blasi & Colafrancesco (1999), in the context 
of the secondary electron model. Lieu et al. (1999) also correctly pointed out that the 
assumption of equipartition is limited by the production of gamma rays through neutral 
pion decay, but this flux was claimed to be much smaller than the EGRET sensitivity, 
falling below the EGRET upper limit already imposed on the gamma ray flux from the 
Coma and Virgo clusters (Sreekumar et al. 1996). 
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We calculate here the flux of gamma rays from the Coma and Virgo clusters in the 
assumption of equipartition of CRs with the thermal energy in the cluster, for two different 
models of the CR injection in the intracluster medium (ICM), and for different injection 
spectra and find that in some cases the gamma ray flux is in excess of the EGRET limit. 
Moreover, we find that for injection CR spectra flatter than E~ 2A (for E 3> m p c 2 ) some 
currently operating experiments like STACEE, HEGRA and Whipple could detect the 
gamma ray signal from Coma and Virgo in the TeV range, provided the CRs are in 
equipartition, or put strong constraints on this condition if no signal is detected. 

The paper is planned as follows: in section ^ we outline the calculations of the gamma 
ray fluxes from clusters; in section |3| we describe the models of CR propagation that we 
used and in section [| we describe our results for the Coma and Virgo clusters. 

2. The gamma ray fluxes 

In this section we calculate the flux of gamma rays due to the decay of neutral pions 
produced in CR collisions in the ICM. This channel provides the dominant contribution to 
gamma rays above 100 MeV. 

Independent of the sources that provide the CRs in clusters, the equilibrium CR 
distribution is some function n p (E p ,r) of the proton energy E p and of the position in the 
cluster. For simplicity we assume the cluster to be spherically symmetric, so that the 
distance r from the center is the only space coordinate. We determine n p for different 
injection models in the next section. The rate of production of gamma rays with energy 
Ery per unit volume at distance r from the cluster center is given by Blasi & Colafrancesco 
(1999) 

g 7 (£ 7 ,r) = 2n H (r)c P dE„ [ ' dEpE^o (E w , E p ) ™ p ^ (1) 
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where E n is the pion energy, E™ m = E 1 + m\j{AE^) is the minimum pion energy needed to 
generate a gamma ray photon with energy E 1 and E^ ax is some maximum energy in the 
injected CR spectrum (our calculations do not depend on the value of E™ ax ). Here uh{t) 
is the density of thermal gas at distance r from the cluster center. For Coma we model the 
gas density through a King profile: 



n H {r) = n 



1 + 



-3/3/2 



(2) 



where ph 400 kpc is the size of the cluster core, no ~ 3 x 10 -3 cm~ 3 and (3 is a 
phenomenological parameter in the range 0.7 — 1.1 (Sarazin 1988) (we use (3 = 0.75). 

For Virgo, we fit the gas density profile given by Nulsen & Bohringer (1995) to find 

n H (r) = 0.076 • ( cm" 3 . (3) 

The function F n in eq. (|l|) represents the cross section for the production of neutral pions 
with energy E n in a CR collision at energy E p in the laboratory frame. Determining this 
function is complicate in the low energy regime where data is scarse. A possible approach 
was proposed by Dermer (1986) and recently reviewed by Moskalenko & Strong (1998) 
and is based on the isobar model. This approach is valid for CR collisions at E p < 3 GeV 
and consists in treating the pion production as a process mediated by the generation and 
decay of the resonance A(1232) in the pp interaction. We refer to the papers by Dermer 
(1986) and Moskalenko & Strong (1998) for the detailed expressions for F w . For E p > 7 
GeV the scaling approach is an excellent approximation of the function F n . In this regime 
the differential cross section for pp collisions can be written as 

(Ia .{E p ,E K ) = ±-a Q f ie {x) (4) 



dE-x En 

where x = E^/Ep, a = 3.2 x 10~ 26 cm 2 and f n (x) = 0.67(1 - x) 3 - 5 + 0.5e~ 18x is the so 
called scaling function. In the scaling regime, the function F w coincides with the differential 
cross section given in eq. (01). 
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Once the gamma ray emissivity is known from eq. ([!]), the flux of gamma rays with 
energy E~ is simply given by volume integration 



^{E 7 ) = -^—p J q dr4vrr 2 g 7 (£ 7 ,r) 



(5) 



where d is the distance to the cluster and R c i is the cluster radius. In fact R c i here plays 
the role of the size of the region where the non thermal processes are observed. We adopt 
here the value suggested from radio observations in Coma, R c i ~ 1 Mpc. This is however a 
very conservative case and it seems likely that magnetic fields extend to larger regions. In 
fact in Ensslin, Wang, Nath & Biermann (1998a) the injection of energy due to formation 
of black holes in the Coma cluster was estimated and compared with the thermal energy in 
a region of Mpc (h^o = h/0.5). For the Coma cluster we shall also consider this less 
conservative case. 



Several sources of CRs in clusters of galaxies were discussed by Berezinsky, Blasi & 
Ptuskin (1997) and it was argued that the known sources (AGNs, radiogalaxies, accretion 
shocks) are not able to provide CRs in equipartition with the thermal gas. An intense and 
short period of powerful emission from the cluster sources was also considered, consistently 
with the observed iron abundance in the cluster, with the same conclusion. Since recent 
observations of non thermal radiation from clusters seem to suggest that equipartition is 
indeed required, we do not make here any assumption on the type of sources and instead 
we assume equipartition and analyze the observational consequences of this assumption. 

The equipartition energy can be easily estimated from the total thermal energy of the 
gas, assuming it has a temperature T: 



3. The Cosmic Ray Distribution 




(6) 
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where njj{r) is given by eq. (Q) for the Coma cluster and by eq. (|) for the Virgo 
cluster. The temperature adopted for Coma is T Co ma = 8.21 k while for Virgo we used 
Tyirgo = 1-8 k (Nulsen & Bohringer 1995). Therefore, from the previous equation we obtain 
E eq = 1.6 x 10 63 erg for Coma and E eq = 1.5 x 10 62 erg for Virgo. These numbers could 
underestimate the total thermal energy due to the contribution of gas out of the ~ 1 Mpc 
region. In fact Ensslin et al. (1998a) estimated for the Coma cluster that the thermal 
energy in a region of 5h§Q Mpc is ~ 1.3 x erg, a factor ~ 6 larger than estimated 

above. They also calculate the expected injection of total (thermal plus non thermal) 
energy due to black hole formation in the cluster, and find in the same region a similar 
number. 

Since not only the energy budget in CRs is not known, but also their spatial distribution 
is very poorly constrained, we consider here two extreme scenarios for the injection of CRs 
and we calculate the equilibrium CR distribution from the transport equation. 

i) Point source 

As argued by Berezinsky, Blasi & Ptuskin (1997), Colafrancesco & Blasi (1998) and 
Ensslin et al. (1998) it is likely that for most of the cluster's age the main contributors to 
CRs in clusters are located in the cluster core. This is the case if a radiogalaxy or more 
generally a powerful active galaxy or a shock produced by merging is the source/ accelerator 
of CRs. There is an additional argument that plays in favor of a source mainly concentrated 
in the center of the cluster: if the average spatial distribution of the galaxies in a cluster is 
not a strong function of time, then it is reasonable to assume that at all times, as today, the 
distribution of the sources is peaked around the cluster center. According with Ensslin et 
al. (1998a) (see also references therein) the spatial distribution of galaxies in Coma is well 
represented by a King-like profile n ga i(r) = [1 + (r/r g ) 2 ]~ - 8 , with r g ~ 160 kpc, appreciably 
smaller that the cluster core, so that a source concentrated in the center seems a reasonable 
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assumption. 

Therefore we assume that the source can be modelled as a point source with a rate of 



injection of CRs given by a power law in momentum Q(E P ) = QoP p ^, where p p = yE 2 — m 2 p 
is the CR momentum and the normalization constant is determined by energy integration 

Qo J Q P dTpTpp' 1 = L p , (7) 

where T p is the kinetic energy and L p is the CR luminosity at injection, forced here to be 
correspondent to the establishment of equipartition in the cluster. We estimate it averaging 
the equipartition energy on the age of the cluster: L p m E eq /to- 

The transport equation that gives the distribution of CRs at distance r from the source 
and after a time t, namely n p (E p , r, t), can be written in the form 

^ p( ^ ,r,t) - D(E p )V 2 n p (E p ,r,t) - [ b ( Ep )n p (E p ,r,t)} = Q(E p )5(r), (8) 

where D(E P ) is the diffusion coefficient and b(E p ) is the rate of energy losses. 

As shown in (Berezinsky, Blasi & Ptuskin 1997, Colafrancesco & Blasi 1998 , Blasi & 
Colafrancesco 1999) for CR protons the energy losses can be neglected and eq. (§) has the 
simple solution (Blasi & Colafrancesco 1999) 

where r max (E p ) = [4D(E p )t] 1 ^ 2 is the maximum distance that on average particles with 
energy E p could diffuse away from the source in the time t. We are interested here in 
the case t = to (to here is the age of the cluster, taken as comparable with the age of the 
universe). The solution of the equation r max (E p ) = R c i gives an estimate of the maximum 
energy E max for which CRs can be considered confined in the cluster volume for all the age 
of the cluster. For reasonable choices of the diffusion coefficient the confined CRs provide 
the main contribution to the energy budget of CRs in clusters and also to the integral 
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flux of gamma rays above 100 MeV, calculated as explained in the previous section. As 
far as gamma rays produced by interactions of confined CRs are concerned the flux of 
gamma radiation is independent on the choice of the diffusion coefficient, as pointed out by 
Berezinsky, Blasi & Ptuskin (1997), and the spectrum of gamma rays simply reflects the 
spectrum of the parent protons (for E 1 > 1 GeV) . Rigorously this is true only for spatially 
constant intracluster gas density, while a density profile, as assumed here, results in a weak 
dependence of the gamma ray spectrum on the diffusion details. Therefore, for the sake of 
completeness we adopt here a specific choice of the diffusion coefficient: we assume that 
the fluctuations in the magnetic field in the cluster are well represented by a Kolmogorov 
power spectrum P(k) oc k~ 5 ^ 3 and we calculate the diffusion coefficient according with the 
procedure outlined by Colafrancesco & Blasi (1998), which gives 

2/3 



D(E P ) = 2.3 x Uf»E p (GeV)V>B?'* 



cm 2 /s (10) 



where is the value of the magnetic field in \xG and l c is the scale of the largest eddy in 
the power spectrum of the magnetic field. 



Eqs. flrO) and (0) completely define the distribution of cosmic rays in the cluster in the 



case of a point source. 

i) Spatially homogeneous injection 

As pointed out above, the budget of CRs in clusters is largely dominated by confined 
CRs, so that in the case of spatially homogeneous injection the distribution of CRs can be 
easily written in the form 



n p (E p ,r) = no-^p p 7 (11) 



ttot 

V 

where V = (4/3)7nR^ is the injection volume, e tot is the total energy injected in the cluster 
in the form of CRs and n is calculated by the normalization condition 

rE rnax 

10 



n 



/ P dT p T pPp "y = E eq , (12) 

J 
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where E eq is calculated according to eq. @. Clearly eq. ([H]) does not describe well the 
CR distribution very close to the cluster boundary. Moreover at sufficiently high energy, 
where CRs are not confined in the cluster volume the CR spectrum suffers a steepening to 
E~^ 1+r, \ with n — 1/3 for a Kolmogorov spectrum. 

4. Results and conclusions 

We study the observational consequences of the assumption of equipartition between 
CRs and thermal gas in clusters of galaxies. In particular we calculated the flux of 
gamma radiation from the Coma and Virgo clusters when equipartition is assumed. This 
assumption seems to be required if a ICS origin is accepted for the hard and soft X-ray 
excess and for the EUV flux from Coma and other clusters of galaxies (note however that 
alternative possibilities can be proposed). In particular, according to Lieu et al. (1999), in 
order to account for the observed cluster soft excess flux from Coma, equipartition between 
CRs and gas is unavoidable. In (Berezinsky et al. (1997)) different possible models of CR 
injection in clusters were considered, including active galaxies, accretion shocks during the 
formation of the cluster and a possible bright phase in the past of the cluster galaxies, 
but none of these sources could account for CR energy densities larger that 1 — 5% of the 
equipartition value, if a conversion efficiency of ~ 10% was assumed for the injection of non 
thermal energy from the total energy of the sources. 

On the other hand Ensslin et al. (1998a) compared the thermal energy in a 5/igo 1 Mpc 
region with the energy injected during the formation of massive black holes in the Coma 
cluster. The total energy (thermal plus non thermal) released in this process was estimated 
to be comparable with the thermal energy in the cluster (if an efficiency factor is assumed, 
the non thermal energy may be smaller than the equipartition value). 
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Since our knowledge of the sources of CRs in clusters is still very poor, we decided to 
adopt here a phenomenological approach and try to find observational tests or consequences 
of our assumptions. The most striking consequence of a large abundance of CRs in a cluster 
is the production of gamma rays through the generation and decay of neutral pions in pp 
interactions. Since the EGRET instrument put an upper limit on the flux of gamma radiation 
above 100 MeV from Coma and Virgo (F EGRET (> lOOMeV) 4 x 10~ 8 phot/(cm 2 s) 
(Sreekumar et al. 1996)), we can use this constraint to test the equipartition assumption. 

Our calculations were carried out for two extreme models of injection of CRs in the 
cluster, namely a point source in the cluster core and a spatially homogeneous injection in 
the cluster volume. In the case of a point source, we can think of it as an effective source, 
in the sense that on average a dominant source or a set of sources are located at the cluster 
center. In this sense it is not needed that the same source remains active for all the age of 
the cluster. 

The energy spectrum of the injected CRs was assumed to be a power law in momentum, 
as expected for a shock acceleration spectrum, and two extreme values of the power index 
were studied, namely 7 = 2.1 and 7 = 2.4, which encompass the whole range of power 
indexes expected from shock acceleration (other models of acceleration also give power laws 
in the same range of parameters). 

Since the gamma ray spectra depend (although very weakly) on the choice of the 
diffusion coefficient, we made here a specific choice, modelling the spectrum of fluctuations 
of the field by a Kolmogorov spectrum and calculating the diffusion coefficient according 
with eq. (|10|). In the numerical calculations we used = 0.1 and l c = 20 kpc (if for 
instance we use = 1 the results on the integral fluxes change only by ~ 10%, confirming 
the weak dependence on the diffusion details mainly due to the use of a specific gas density 
profile). 
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The integral fluxes of gamma radiation above 100 MeV for the cases mentioned above 
and in the conservative scenario of R c i = 1 Mpc, are reported in Table p] for the Coma and 
Virgo clusters. 

EDITOR: PLACE TABLE g HERE. 

Due to the appearence of the flat region at low gamma ray energy, typical of spectra 
from pion decay, there is not a strong dependence of the integral flux on 7. In some cases 
considered the gamma ray flux exceeds the EGRET upper limit. As it could be expected, 
the gamma ray flux is larger for the case of a point source in the cluster center and the 
EGRET limit is exceeded by a factor ~ 1.7 for Coma and by a factor ~ 9 for Virgo. In the 
case of homogeneous injection the gamma ray fluxes are slightly smaller than the EGRET 
upper limit both for Coma and Virgo. 

These results are more impressive when the condition of equipartition is imposed on a 
larger region of size Mpc (Ensslin et al. (1998a)): for a single source in the cluster 
center the EGRET limit is exceeded by ~ 9 for 7 = 2.1 and by ~ 8 for 7 = 2.4. For an 
homogeneous injection the predicted fluxes are in excess of the EGRET limit by ~ 7 for 
7 = 2.1 and by ~ 6 for 7 = 2.4. 

It is worthwhile to stress again that this result is practically independent on the 
specific choice of the diffusion coefficient. In fact, the CRs relevant for the production of 
gamma rays in the energy range 0.1 — 10 GeV are certainly confined in the cluster for any 
reasonable choice of the diffusion coefficient, and the spectrum in this region is independent 
on this choice. The effects of the diffusion may appear only at higher energy where gamma 
rays are produced by CRs not confined in the cluster. As shown in (Berezinsky et al. 1997) 
this produces a steepening of the gamma ray spectrum to a power law with index 7 + 77 
(with 1] = 1/3 in the case of a Kolmogorov spectrum of fluctuations) at energies larger than 
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a knee energy Ex- At smaller energies the gamma ray spectrum reproduces the spectrum of 
the parent CRs. The transition appears at Ex oc R 2 c J (B^'Hq) , as obtained from the 
equation r max (E p ) = R c i. Actually this was shown in (Berezinsky et al. 1997) for the case 
of a constant intracluster gas density. In the more realistic case considered here, where the 
gas is modelled by a King or a power law profile, the gamma ray spectrum suffers a smooth 
steepening even for confined CRs, but this affects the integral flux above 100 MeV only at 
the level of ~ 10%. 

The integral spectra of gamma rays from Coma with energy > E 1 as functions of the 
energy E 1 are shown in Fig. la (for the point source) and lb (for the homogeneous case) 
for R c i = 1 Mpc. In the same plot we draw the sensitivity limits for several present and 
planned experiments for gamma ray astronomy. The solid lines refer to 7 = 2.1 while the 
dashed lines are obtained for 7 = 2.4. 



EDITOR: PLACE FIGURE [l] HERE. 



The fluxes in the energy region E 1 < 100 GeV are well above the detectability limit 
of GLAST, so that there is no doubt that the question of equipartition will be completely 
answered with the next generation gamma ray satellites. However fig. 1 also shows that the 
signal from Coma could be detectable even in some current experiments, provided 7 < 2.4. 
In particular STACEE could detect the signal above 30 GeV and Whipple might detect 
the signal for E 1 > 250 GeV. The flux should be detectable by the HEGRA Cerenkov 
telescope above 500 GeV. A non-detection from these experiments would imply a reduction 
of the energy density in clusters by about one order of magnitude below equipartition for 
7 = 2.1. For steep spectra only STACEE has a slim chance to detect the signal. In the 
same energy range the next generation gamma ray experiments will very likely measure the 
flux of gamma rays for any value of 7 in the range considered here. 
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In the case of the Virgo cluster and R d = 1 Mpc, the fluxes are plotted in fig. 2a (for 
sources in the center) and fig. 2b (for a homogeneous injection) and conclusions similar to 
the ones outlined for Coma hold. 

Note that this result is subtantially different from the one obtained in previous 
calculations. In particular Ensslin et al. (1997) reached the conclusion that the fluxes from 
Coma and Virgo are orders of magnitude too low to be detectable in the TeV range. This 
conclusion was obtained because following Dar and Shaviv (1995) the gamma ray spectrum 
was assumed to reproduce the equilibrium spectrum of CRs in the Galaxy oc E~ 2 ' 7 (this did 
not affect appreciably their integral fluxes above 100 MeV, which are not very different from 
the ones obtained here for the homogeneous case). However, as shown in (Berezinsky et al. 
1997, Colafrancesco & Blasi 1998, Blasi & Colafrancesco 1999 ) and confirmed here, the 
spectrum of gamma rays from pp collisions in clusters does not reproduce the equilibrium 
CR spectrum , but the generation spectrum as far as gamma ray photons are produced by 
interactions of CRs confined in the cluster, as it is the case for gamma rays with energy 
less than ~ 1 — 10 TeV, for the values of the parameters used here (though, as pointed out 
before, a slight steepening is introduced by the gas density profile). Therefore the gamma 
ray spectrum from CRs in clusters is approximately E~^ up to some maximum energy Ek 
where CRs begin to be no longer confined in the cluster volume. As a consequence the 
gamma ray fluxes in the TeV range could be detectable even by present experiments if 
the CRs are in equipartition with the gas in the cluster. On the other hand, if no flux is 
detected, this will put a strong constraint on the equipartition assumption. 

While the low energy integral gamma ray flux is very weakly dependent on the choice 
of the diffusion coefficient, the correspondent flux at higher energies is more sensitive to 
it, since, as explained above, the position of the knee is affected by this choice. In the 
context of a Kolmogorov spectrum, the maximum diffusion coefficient is obtained for a 
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larger value of l c . The choice l c m 20 kpc was inspired by the typical size of the galaxies in 
the cluster. The largest scale where the magnetic fluctuations are injected is the typical 
distance between galaxies, of the order of l c ~ 100 kpc. For l c « 20 kpc the position of the 
knee is at ~ 1 — 10 TeV, while for l c 100 kpc the knee is at ~ 10 — 20 GeV. However, 
since the steepening in the gamma ray spectrum begins at large energy, the difference in the 
plots caused by the use of this larger diffusion coefficient is a factor ~ 2 at ~ 100 GeV (for 
R c i = 1 Mpc), so that the possibility of detecting the gamma ray fluxes in this energy region 
is not appreciably affected and remains an interesting possibility. The situation improves 
rapidly with an increasing R c i. In fact the value of the knee energy for a Kolmogorov 
spectrum goes like R® cl and high energy CRs are easily confined in a region of 4 — 5 Mpc. If, 
following Ensslin et al. (1998a), we use R c i = 5h^ Mpc, then, as shown before the absolute 
gamma ray fluxes increase at all energies by a factor ~ 5 — 10 and the steepening at high 
energy is only found at E ^ 10 3 TeV for the diffusion coefficient in eq. (JTU ) . 



At present we can only use the EGRET limit as a constraint. For Coma this limit 
implies that the energy density must be smaller than ~ 60% of the equipartition value if 
the CRs are mainly contributed by sources in the central part of the cluster (with R c i = 1 
Mpc). If injection occurs uniformly over the cluster volume, than the equipartition CR 
energy density is compatible with the EGRET limit (for R c i = 1 Mpc). For Virgo cluster 
the EGRET limit implies that the energy density in CRs must be smaller than ~ 10% of 
the equipartition value for the case of sources in the center. 

As stressed above, the case of equipartition in a larger region, (Ensslin et al. (1998a)) 
is already ruled out by present gamma ray observations by EGRET. As a consequence the 
CR energy density in this case for Coma is forced to be ^ 12% of the equipartition value 
for a point source and ^ 14% for homogeneous injection. 



EDITOR: PLACE FIGURE @ HERE. 
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We suggest that experiments like STACEE, HEGRA, Whipple and future gamma ray 
experiments look at the signal from nearby clusters, because this could definitely confirm or 
rule out the possibility that equipartition of CRs with the thermal gas is achieved in clusters 
of galaxies, or at least impose new and stronger constraints on the maximum allowed CR 
energy density in clusters. 

The author is grateful to A. Olinto, S. Colafrancesco, C. Covault and R. Ong for many 
useful discussions and to the anonymous referee for several interesting comments. The 
research of P.B. is funded through a INFN fellowship at the University of Chicago. 
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Table 1: Summary of the gamma ray integral fluxes at E 1 > 100 MeV 



Model 


<$>Coma{E > lOOMeV) phot/{cm 2 s) 


$ V irgo(E > lOOMeV) phot/(cm 2 s) 


7 = 2.1, Point source 


7 x 10~ 8 


3.7 x 10~ 7 


7 = 2.4, Point source 


6.5 x 10~ 8 


3.6 x 10~ 7 


7 = 2.1, Homogeneous 


2.5 x 10~ 8 


3 x 10~ 8 


7 = 2.4, Homogeneous 


2.1 x 10~ 8 


2.6 x 10' 8 



-20- 



Fig. 1. — Integral fluxes of gamma radiation from the Coma cluster compared with the 
sensitivity limits of some current and proposed gamma ray experiments (thick solid lines). 
The thin solid line curves are referred to 7 = 2.1 and the dashed ones to 7 = 2.4. In both 
cases R c i = 1 Mpc. The EGRET limit at 100 MeV is indicated by the arrow, a) Injection 
of CRs by a point source in the cluster center, b) Homogeneous injection of CRs over the 
cluster volume. 

Fig. 2. — The same as Fig. 1 for the case of the Virgo cluster. 



